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ABSTRACT 

Aims. We want to study the IR (> 8 yum) emission of samples of galaxies selected in their rest-frame UV in a very homogeneous 
way (wavelength and luminosity) from z = to z = 1 : comparing their UV and IR rest-frame emissions will allow us to study the 
evolution of dust attenuation with z as well as to check if a UV selection is able to track all the star formation. This UV selection will 
also be compared to a sample of Lyman Break Galaxies selected at z — 1 • 

Methods. We select galaxies in UV (1500-1800 A) rest-frame at z = 0, z = 0.6-0.8 and z = 0.8-1.2 together with a sample of 
Lyman Break Galaxies at z = 0.9-1.3, the samples are built in order to sample the same range of luminosity at any redshift. The UV 
rest-frame data come from GALEX for z < 1 and the f/-band of the EIS survey (at z = 1). The UV data are combined to the IRAS 
60 yum observations at z = and the SPITZER data at 24 fan for z > sources. The evolution of the IR and UV luminosities with z is 
analysed for individual galaxies as well as in terms of luminosity functions. 

Results. The L IR /L UV ratio is used to measure dust attenuation. This ratio does not exhibit a strong evolution with z for the bulk of 
our sample galaxies but some trends are found for galaxies with a strong dust attenuation and for UV luminous sources: galaxies with 
L IR /L UV > 10 are more frequent at z > than at z = and the largest values of L 1R /L V \j are found for UV faint objects; conversely the 
most luminous galaxies of our samples (L uv > 2 x 10'°L o ), detected at z = 1, exhibit a lower dust attenuation than the fainter ones. 
L IR /L UV increases with the K rest-frame luminosity of the galaxies at all the redshifts considered and shows a residual anti-correlation 
with L uv . The most massive and UV luminous galaxies exhibit quite large specific star formation rates. Lyman Break Galaxies exhibit 
systematically lower dust attenuation than UV selected galaxies of same luminosity but similar specific star formation rates. 
The analysis of the UV+IR luminosity functions leads to the conclusion that up to z = 1 most of the star formation activity of UV 
selected galaxies is emitted in IR. Whereas we are able to retrieve all the star formation from our UV selection at z = 0.7, at z = 1 we 
miss a large fraction of galaxies more luminous than ^ 10 U L Q . The effect is found larger for Lyman Break Galaxies. 

Key words, galaxies: evolution-galaxies: stellar content-infrared: galaxies-ultraviolet: galaxies 



1. Introduction 

The measure of the star formation rate at various redshifts is 
one of the most spread out diagnostic to quantify the evolution 
of the galaxies, individually or for the population as a whole. 
The best way to perform such an analysis is to rely on galaxy 
samples selected in the same way at different redshift. As long 
as star forming galaxies have to be selected a UV selection is 
theoretically very efficient. Nevertheless, the situation is made 
difficult because of dust attenuation. Indeed dust attenuation 
affects the measure of the star formation rate derived from the 
observed UV emission and a correction (sometimes quite large) 
must be applied to the observed UV emission before translating 
it into SFR. However the effects of dust obscuration may be 
even more dramatic if they lead to a loss of galaxies which 
might not be detected at all in UV surveys: in such a case a 
correction for dust attenuation of the total light observed in UV 
would not be sufficient to retrieve the whole star formation at a 
given redshift. 



The UV (1500-1800 A) to IR (8-1000 yum) luminosity 
ratio L IR /L UV is now commonly used as a robust proxy for dust 
attenuation. The GALEX all sky survey associated with the 
IRAS catalogues has produced large samples of nearby galaxies 
observed at both wavelengths and used to study the variation of 
Lir/Luv with the total Lir + Luv luminosity (i.e. the total SFR 
of galaxies). Lir/Luv has been found to increase with Lir + Luv 
(Mar tin et aT| |2005t[Buat et al.||2007| l. The same relation hold s 



for galaxies selected either in UV or in IR (Buat et al. 2007). 



This confirms the general increase of dust attenuation with the 



luminosity of the galaxies already reported (e.g. Hopkins et al. 
[200T]|Moustakas et al.||2006] i. 

Such a relation between total luminosity and dust attenuation 
in galaxies is useful to correct for systematic effects of dust 
attenuation in galaxy surveys. It is important to check if this 
relation still holds at higher z. Several studies have been devoted 
to this issue. Contrary to what is found in the nearby universe, it 
appears that the result depends on the way galaxies are selected. 
IR selected galaxies at intermediate redshift (z - 0.5-0.8) seem 
to follow the mean trends found at z = between Lir + Luv 
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and Lir/L uv (Choi et al. 2006 Xu et al. 2007, Zheng et al. 
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2006j). |Buat et al.| (|2007) reported only a slight decrease of dust 
attenuation (~ 0.5 mag) for luminous IR galaxies (LIRGs) at 
Z = 0.7 as compared to a similar sample of galaxies at z = 0. 
When galaxies are selected in UV/optical the situation is quite 
different: it has been reported a strong decrease of Lir/Luv for 
a given total Lir + Luv as compared to what is found at z = 
HBurgarella etaEl [200r3l [20071 |Reddy et al.||2008| l. 
The origin of this discrepancy is not clear. Are the properties 
of galaxies changing with z and/or do we sample very different 
galaxy populations when selecting in UV or in IR? In the 
nearby universe |Buat et aL ( |2007] > have shown that a UV or 
an IR selection leads to similar results except for intrinsically 
very luminous galaxies which are under-represented in the UV 
samples. Such objects are rare at z — but because of the 
evolution of the luminosity functions with z we might expect a 
different situation at higher z. 

In order to answer to these questions, we gather here several 
samples of UV selected galaxies from z — to z — 1.2 selected 
in a very homogeneous way. We also consider a sample of 
Lyman Break Galaxies (LBGs) selected at z — 1 to be compared 
to our pure UV selections. We add IR fluxes (when available) 
for all the galaxies of these samples. With these data we want 
to study the evolution with z of the IR emission mainly through 
the study of the Lir/Luv ratio and to measure the total star 
formation activity by combining IR and UV emissions. The 
reliability of UV selected galaxies to trace all the star formation 
will be discussed through the analysis of bolometric (Lir + Luv) 
luminosity functions. 

Throughout the paper we will assume Q,„ = 0.27, Qa = 0.73 
and Hq = 71 km s _1 Mpc~'. All the magnitudes are given in 
the AB system except for the R magnitude from the COMBO- 
17 survey (section 2.1). The luminosities are defined as vL v and 
expressed in solar units assuming L = 3.83 x 10 33 erg s" 1 



2. The galaxy samples 

2.1. The UV selected samples 

The samples must be purely UV selected. As a consequence we 
mostly rely on the GALEX survey for redshift lower than 1 . 
Atz = 0, we take the IRAS/GALEX sample built by |Buat et aT] 
(2007 ). This sample consists of galaxies selected in the GALEX 
FUV band (1530 A) and with FUV < 17mag. It is a flux lim- 
ited sample and the luminosity function has been built down to 
Lfuv = 1O 8 L . The FUV wavelength will be taken as the refer- 
ence wavelength for the samples at higher z, and quoted as UV 
throughout the paper. 

At higher z, our samples are extracted from the GALEX 

2005) 



deep observations of the CDFS. GALEX (jMorrissey et al. 



observed this field in both the FUV (1530 A) and the NUV 
(2310 A) as part of its deep imaging survey. In order to add IR 
data to the galaxy sample we restrict our study to the sub-field 
covered by SPITZE R/MIPS observa tions as part of the GOODS 
key program (e.g. Elbaz et al. 2007 1. 



Th e sample at z — 0.7 was already used by Buat et al. 
\200S\ , at this redshift the NUV band of GALEX at 2310 A 
corresponds to the FUV rest frame of the galaxies. This sam- 
ple is thus made of galaxies selected in NUV and with a red- 
shift comprised between 0.6 and 0.8. Redshifts come from the 
COMBO-17 survey of the field ( |Wolf et al.j [2004| >. The reduc- 
tion of the data and the cross-identification with the COMBO- 



son between the fluxes given by the GALEX pipeline, aper- 
ture photometry measurements and fluxes obtained by PSF fit- 
ting with DAOPHOT. From this comparison we conclude that 
the fluxes given by DAOPHOT and previously used might be 
slightly over-estimated by 0.2 mag. Although it is not clear if 
this shift is reliable for crowded objects we decided to apply the 
correction to all the UV data in the CDFS field: NUV(new) = 
NUV(old) - 0.2 mag. This systematic correction remains very 
small (of the order of the error) and does not imply any modifi- 
cation of our previous results. Whereas the completeness of the 
NUV data at a level of 80 % is obtained for NUV = 26 mag, 
we truncate the sample at NUV = 25.3 mag in order that more 
than 80 % of the GALEX sources are identified in COMBO-17 
with R < 24 mag. This limit ensures a redshift accuracy better 
than 10% (Wolf et al.]|2004) . We also restrict the final sample to 
objects with a single counterpart in COMBO-17 within 2 arcsec 
(i.e. 90 % of the UV sources). 300 galaxies are thus selected. The 
limit of NUV = 25.3 mag corresponds to log Luv = 9.3(L Q ) at 
z = 0.7. 44 % (131/300) of these sources are detected at 24 fim. 
For the undetected ones, we adopt an upper limit at 0.025 mJy 
dBuatet al.||2008l >. 

At z — 1 we must build a new sample. The NUV band of 
GALEX corresponds to 1 155 A in the rest frame of the galaxies. 
Downward 1200 A the spectral energy distribution of galaxies 
is poorly known and the few available observations have shown 
that the shape of the spectral energy distribution (SED) may vary 
a lot from galaxy to galaxy (Bu at et al.| |2002| |Leitherer et al.| 
|2002 ): the power-law model valid for A > 1200 A cannot be 
safely used at shorter wavelengths. To avoid this difficulty we 
perform a U selection. We again work in the CDFS/GOODS 
field and f/-band observations have been performed as part of 
the EIS survey ( Arnouts et al. 200 1 1 > . We have cross-correlated 
the U selected catalogue with the COMBO-17 sample with a 
tolerance radius of 1 arcsec. Redshifts are from COMBO- 17+4 
(Tapken, private communication) that provide safe redshifts, es- 
pecially at z > 1 in the GOODS-South field by combining 
COMBO-17 filters with three near infrared bands in ISAAC JHK 
bands. 69 % of the U sources are identified uniquely down to 
U = 26 mag. We cut the sample at U — 24.3 mag which corre- 
sponds to 80 % of the sources identified in COMBO-17 with a 
single counterpart. The photometric redshifts from COMBO-17 
are secure only for galaxies with R < 24 mag (in the Johnson 
system). 96 % of our galaxies fill this condition. We then se- 
lect galaxies in the redshift bin 0.8-1.2 which results in a sam- 
ple of 316 galaxies, at this redshift the L-band corresponds to 
~ 1800 A in the galaxy rest-frame (z = 1) when our refer- 
ence wavelength is 1530 A. Because of the uncertainties about 
the shape of the UV SEDs we prefer to avoid any interpolation 
between the NUV and U observed fluxes and keep the uncor- 
rected U data. We can estimate the uncertainty due to the shift 
of the rest-frame wavelengths from 1800 to 1530 A by assuming 
a power-law for the continuum between these two wavelength, 
for reasonable values of the power-law valid for a UV selection 
(a — -2 to 1 with F a oc A a ) the error is at most 20 %. The limit- 
ing magnitude U = 24.3 mag corresponds to log Luv = 9.9 (L Q ) 
at z — 1. As for the sample at z — 0.6-0.8, the 24 /im data from 
SPITZER/MIPS were cross-correlated to the U sources identi- 
fied in COMBO-17 within a tolerance radius of 2 arcsec dBuatet] 



17 sources is also described in |Burgarella et al. (2006i. For the 
brightest isolated sources of the field we perform a compari- 



|aL)|2008| ). 207 out of the 316 galaxies are detected at 24 fim. For 
the undetected ones, we again adopt an upper limit at 0.025 mJy 
dBuatet al.||2008l >. 

The aim of this paper is to determine if a UV selection is able 
to track all the star formation. The comparison of the IR and UV 



V. Buat et al.: UV selected galaxies from z = to 1 



3 



emissions will allow us to study star formation and dust atten- 
uation affecting newly formed stars. In this context quasars and 
active galaxies are excluded from all our samples. We exclude 
objects classified as QSO/Seyfertl in the Combo-17 classifica- 
tion (2 objects at z — 0.7 and 4 at z — 1). We have also cross- 
correlated our galaxy sample with the X ray sources observed 
by CHANDRA in the field ( [Bauer et al.l[2004| and we discard 
22 sources at z — 0.7 and 19 at z — 1 as X-ray emitters. As a 
final check we can compare mid-IR IRAC colours as suggested 
by Stern et al. ( 2005 | l to determine if our samples were still con- 
taminated by AGNs. 1 and 4% of the samples at z=0.7 and 1 
respectively were found in the AGN area close to the boundary 
where the contamination by star forming sources is significant. 
We choose not to exclude these remaining sources. 



2.2. Lyman Break Galaxies atz-\ 

This paper presents a sample of galaxies in the CDFS-GOODS 
field selected as being LBGs with a similar UV selection as that 
used in previo us works (fB urgarell a'et al.| ( |200"6] ), [Burgarella et al. 
( 2007 1 and Burgarella et al.| ( |2009| l). The primary selection was 
performed in GALEX NUV band down to NUV = 25.3 mag, 
as discussed in Section 2.1 this limit ensures us that 80 % of 
the sources are identified with a reliable redshift in COMBO-17. 
Then we look up into the GALEX FUV band for counterparts 
down to FUV = 26.8 mag (80 % completeness on the GALEX 
detections) . To be selected as a Lyman Break Galaxy, an object 
must comply with the two criteria: 1) its redshift must be in the 
range 0.9 < z < 1.3 and 2) its UV colour FUV - NUV > 2 mag. 
Indeed, some NUV faint objects that would be classified as 
LBGs if the FUV limiting magnitudes objects had reached down 
to FUV < 27.3 mag are not selected as LBGs in the present 
sample. 

This selection provides 1 17 LBGs with a unique counterpart 
in the optical out of which 58 are detected at 24 /mi, for the 
other ones an upper limit of 0.025 mJy is adopted as for the other 
samples. In addition, we have 33 LBGs with two counterparts in 
the optical out of which 20 are detected at 24 /mi. 

As for the z — 1 sample, the UV luminosity is based on the 
observed f/-band. At z = 1 the selection of galaxies on their 
Lyman break focuses on galaxies with a high intrinsic UV con- 
tinuum and a large intrinsic break since the role of the intergalac- 
tic medium on the amplitude of the break is known to be low at 
this redshift (e.g. Malkan et al. 2003[ ). The situation is different 
for LBGs selected at higher z whose break is dominated by the 
effect of the intergalactic medium. Nevertheless Burgarella et al. 
( 2007| l found similar spectral energy disributions for LBGs at 
z - \ and z - 3. 



observations of nearby galaxies by IRAS and/or ISO (e.g.|Chary| 
|& Elbaz| [200T] |Dale & Heloul |2002l [Takeuchi et al.[ |2005a) . 
The recent observations of SPITZER have stressed the large va- 
riety of IR spectral energy distributions in nearby galaxies ( |Dale| 



2.3. Estimating the IR luminosities 

The aim of the paper is to compare the IR and UV emissions of 
the galaxies from z = to z - 1. At z = 0, th e IR (8-1000 /mi) 
emission is estimated using the calibration of |Dale et al.| ((2001 ) 
based on the 60 and 100 fim fluxes from IRAS ( |Buat et al.||2005) . 
At higher z, we must use the emission at 24 /mi to estimate the 
total far infrared emission. At z = 0.7 (resp 1) the observed 
24 /mi corresponds to rest-frame 15 /mi (resp 12 /mi). The 
extrapolation from the mid infrared (MIR) emission to the total 
IR one is known to be quite difficult but up to z — 1 we still 
remain in a wavelength range observed in the nearby universe 
by either IRAS or ISO. As a consequence several calibrations 
and IR SED templates were proposed in the literature based on 



et al. 2005 Rieke et al. 2008) >, these observations imply large 



uncertainties in the extrapolation from a monochromatic flux to 
the total IR emission. To illustrate this uncertainty, we can com- 
pare several calibrations of the monochromatic luminosities at 
12 and 15 /zm into bolometric IR luminosities. We consider the 
calibrations of |Chary & Elbaz| ( |200T] ) and |Dale & Helou|p002 1 
based on templates built by combining IRAS and ISO data on 
small samples of galaxies. The Dale & Helou (2002) templates 
are calibrated in total IR luminosities following the method of 
Marcillac et al. (|2006 )) and we obtain the following relations 



between total and monochromatic luminosities: 
logLi R = 1.25 log L n -0.341 
logL IR = 0.985 logL I5 + 1.26 



(1) 
(2) 



We also consider the calibrations of Tak euchi et al.| (|2005a) 
(hereafter TBI05) based on a statistical analysis of all the galax- 
ies observed in the four bands of IRAS in the Point Source 
Catalog. The relation obtained by TBI05 at A — 12 //m for local 
galaxies (i.e. when we observe a galaxy at /l b s , its emitted wave- 
length A em corresponds to 12 /mi at z = 1) is slightly modified 
to take into account a slight non-linearity for the most luminous 
objects: 

log L m = 2.265 + 0.707 log L n + 0.014(log L n f , (3) 

whereas at A — 15 /mi, we used the original relation 

logL IR = 1.23 + 0.972 log L 15 (4) 

at z = 0.7, since the non linearity is not significant. We have 
checked that the estimated L JR s do not differ by more than 5% as 
compared to the old calibration presented in TBI05. We also add 
to the comparison the recent calibration obtained by |Rieke et al.| 
(2008) at 12 /mi and based on a compilation of SPITZER data 
with a special emphasis on LIRGs and ULIRGs. The compari- 
son between all these calibrations is shown in Fig. [T] In order 
to highlight the differences we have normalized the relations to 
those of DH02. The rms dispersion is overplotted for each rela- 
tion. 

At 15 /mi (corresponding to the sample at z = 0.7 in the present 
wor k) the calibrations of|Chary & Elbazl ( p50T] l (hereafter CE01) 
and [Dale & Helou| ( |2002[ ) (hereafter DH02) are found very sim- 
ilar: the calibration of TBI05 leads to a slightly lower L IR by 
about 0.1 dex than those obtained with the CE01 relation. This 
difference is much smaller than the intrinsic dispersion of the 
correlations and |Buat et aT] ( |2007| l have found that both cali- 
brations led to very similar results about dust attenuation for 
a sample of LIRGs. However, at 12 /mi (z - 1) the discrep- 
ancy is found larger between the calibrations of CE01 and DH02 
on the one hand (again very similar) and that of TBI05 on the 
other hand: it reaches 0.2 dex for galaxies with Li R 10 n L Q 
the TBI05 relation giving systematically lower IR luminosi- 
ties. Nevertheless these relations calibrated at z = remain 
marginally consistent within one rms and we expect at least the 
same amount of dispersion at higher redshift. The relation pro- 
posed by Rieke et al. ( 2008 ) for log Li R > 8.5(LO) appears much 
steeper than the other ones leading to larger IR luminosities es- 
pecially for intrinsically luminous objects. Rieke et al. (2008) 
have gathered SEDs of nearby LIRGs and ULIRGs together with 
the |Dale & Helou| ( [2002| templates applied to the SINGS sample 
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Fig. 1. Different calibrations of log(L24) at z = 0.7 (upper 
panel)- and z — 1 (lower panel) versus log Ltr. All the calibra- 
tions are normalized to the values found with the DH02 relations 
given in the text. The three lines plotted for each calibration cor- 
respond to the mean relation and the rms dispersion. 



for galaxies of intermediate luminosity. Therefore we can infer 
that the discrepancy found between the Dale & Helou ( 2002[ > 
and Rieke etaL] ( |2008| l calibrations is due to the introduction of 
these LIRGs and ULIRGs. Once again these differences illus- 
trate the uncertainty inherent to these calibrations. At z — 0.7, 
26% of our galaxies detected at 24 /itn and 11% of the whole 
sample are LIRGs-ULIRGs; at z 1 these fractions reach 34% 
of the galaxies detected at 24 fim and 22 % of the whole sample. 
Although these fractions are significant, LIRGs-ULIRGs do not 
dominate our sample and we will not use the relation of Rieke 
et al. ( 2008| l. We perform all the analyses reported in this work 
for the two calibrations TBI05 and CE01. The plots are quali- 
tatively similar and they will be presented with the TBI05 cali- 
bration. When quantitative evaluations are made (regressions or 
percentages) they will be given for both calibrations. 



Fig. 2. log(LiR/Luv) versus log(LiR + Luv) for the different sam- 
ples (i.e. redshifts) defined in this work and the BM/BX galaxies 
of Reddy et al. ( 2006| l. The different symbols corresponds to de- 
tections at 24 /im (z = 0: plus, z = 0.6-0.8: circles, z = 0.8-1.2: 
triangles, LBGs: crosses, BM/BX: stars), dots are for upper lim- 
its. The solid line corresponds to the adopted limit in UV lumi- 
nosity: log(L uv [L ] = 9.9). 



3. Variation of Lt R /Ltjv 

As already underlined in the introduction Lir/Luv ls a robust in- 
dicator of dust attenuation as long as we are dealing with galax- 
ies forming stars actively. Reddy et al. ( 2008| > and Burgar eTla et| 
[aL ( 2009) have reported a clear decrease of this ratio for LBGs 
at z = 1 and BM/BX galaxies at z = 2 for a constant La + Luv 
luminosity. Such a decrease of dust attenuation as redshift in- 
creases may have large consequences on the search of high red- 
shift galaxies and the measure of their star formation rate. Here 
we re-investigate this question with our homogeneous samples 
selected in a similar way at z = 0, - 0.7 and 1 . 



3.1. Lir/Ltjv versus L IR + L uv 

As detailed in the previous section the samples are all selected 
in UV rest-frame at wavelengths sufficiently close to avoid K- 
corrections. Nevertheless, the range of luminosities observed in 
each sample is different. Before any comparison we must cut 
the samples at the same luminosity. The most stringent limit 
is for the {/-selected sample. With a cut at U — 24.3 mag 
we are only able to observe galaxies with log L uv > 9.9 (L Q ) 
at z = 1. In Fig. [2] are reported all the galaxy samples con- 
sidered in this work and truncated at log Luv > 9.9 (L Q ). A 
trend with redshift seems to be present with galaxies shifted 
toward the right when z increases. The mean values found for 
Lir/Luv in the luminosity bin 10.8-1 1.3 (L ) are < Lir/Luv >= 
0.86±0.35, 1.07±0.30, 0.75+0.30 and 0.54+0.35 at z = 0, 0.7, 1 
and for the LBG sample respectively. Although no clear trend 
is confirmed from z=0 to 1 given the large dispersion of the 
Lir/Luv distribution, for a given L JR + Luv, dust attenuation as 
traced by Lir/Luv is lower for LBGs at z — 1. 
We have also over plotted the sample of BM/BX galaxies of 
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Fig. 3. log(LiR/Luv) versus log(LiR + Luv), 3 sub-samples are 
defined according to the UV (rest-frame) luminosity and differ- 
ent symbols are used for each of them. The lines correspond to 
the locus of galaxies with a fixed Luv: log Luv = 9.9 (solid line), 
10.1 (dotted line), 10.3 (dashed line) and 10.5 (dot-dashed line) 
in solar units. For the sake of simplicity upper limits at 24 fim 
are not over-plotted 



Fig. 5. fraction of galaxies with log(LiR/Luv) > 0.7, 
log(LiR/Luv) > 1- The symbols are the same as in Fig. [2] The 
first bin of luminosity in the lower panel is not considered at 
z = 0.7 since the upper limits of log(LiR/Luv) calculated for the 
galaxies with this luminosity and undetected at 24 [tm may be 
higher than 0.7. 
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Fig. 6. fraction of LIRGs, the symbols are the same as in Fig. 2] 
The solid lines refer to the TBI05 calibration for Lir and the 
dotted lines to the CE01 calibration 



Fig. 4. log(LiR/Luv) versus Lyy for the different samples de- 
fined in this work. The diagonal solid lines are the limits above 
which galaxies are LIRGs or ULIRGs. 



Reddy et al. (2006) at z=2. The IR luminosities are estimated 



from the 24 fim fluxes using the calibration of Caputi et al 
( |2007| l as in |Reddy efal] ( |2008| l and the UV luminosity is cal- 
culated with the G band fluxes. Given the differences between 
these estimates and those performed for the samples studied in 
this work we do not make a quantitative comparison. Anyway 



the BM/BX galaxies continue the trend reported for LBGs at 

However the interpretation of this plot is made difficult by 
the fact that the quantities reported on the axes are both a com- 
bination of Lir and Luv- This is illustrated in Fig. [3] where the 
galaxies are considered according to their UV luminosity, what- 
ever their redshift is. The lines represent the locus of galax- 
ies with a given Luv- The locus of the galaxies in this plot 
is strongly constrained by their UV luminosity: when Lir/Luv 
varies, galaxies with a given Luv move along lines like those 
over-plotted on the diagram. Therefore the shift seen in Fig. [2] 
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is due to the presence of more luminous galaxies in UV when 
Z increases as expected from the evolution of the UV luminos- 
ity function ( |Arnouts et"aT] |2005| |Takeuchi efaL] |2005b| ). The 
variation of dust attenuation as traced by Lir/Luv can only be 
quantified along these lines of constant Luv- 

3.2. L IR /L UV versus L uv 

From the above analysis, it is clear that we must avoid to com- 
bine IR and UV luminosities on both axes since the resulting 
plot is too strongly constrained. We can analyse the variation of 
Lir/Luv as a function of Luv alone. In such a plot we are not 
affected by volume effects since all the samples are selected in 
UV. In Fig. |4] is reported the variation of Lir/Luv as a func- 
tion of L UV - The loci of LIRGs (L 1R > 10 11 L ) and ULIRGs 
(Lir > 10 12 L Q ) are also indicated. This time the general shape 
of this diagram is strongly constrained by the upper limits at 
24 pm which hamper any discussion about the low values of 
Lir/Luv- The upper envelope of the distribution shows a trend: 
as Luv increases the maximum value of Lir/Luv decreases from 
~ 50 to ~ 3 (~ 65 to ~ 5 for the CE01 calibration). A quan- 
titative interpretation of this varying upper limit is difficult be- 
cause of the statistics: since the total number of galaxies per lu- 
minosity bin decreases as the UV luminosity increases, we ex- 
pect less extreme cases even with a similar parent distribution 
for L IR /Luv- Anyway the most UV luminous galaxies exhibit a 
very moderate dust attenuation (log(LiR/Luv) = 0.5 corresponds 
1 .2 mag with the calibration of |Buat et al. i 



to A 



uv 



( |2005| »and 



aintest ones 



the galaxies with the largest dust attenuation are the 
in UV. 

In order to go further in the interpretation of Fig. [4] we have cal- 
culated the fraction of galaxies with log(L IR /Luv) larger than 0.7 
and 1 (corresponding to Auv = 1.5 and 2 mag respectively (Buat 
|et al.| |2005| l) as well as the fraction of LIRGs for each redshift, 
as a function of Luv- The cuts adopted for log(L IR /Luv) (0.7 
and 1) are chosen to be not affected by the non detections at 24 
fan in the high redshift samples (cf the upper limits reported in 
Fig.|4jl . If the CE01 calibration is used instead of the TBI05 one 
the cuts in log(Li R /Luv) have to be increased of 0.1 dex (i.e., 
log(L IR /Luv) =0.8 and 1.1). 

The results are reported in Fig. [5 and [6] At z = all the 
galaxies are detected in IR and the distribution of Lir/Luv is 
well described by a Gaussian with a mean value of 0.55 dex and 
a standard deviation of 0.3 dex. If we first consider the fraction 
of galaxies with log(LiR/Luv) > 0.7, this fraction is not found to 
be very dependent on the redshift or the UV luminosity at least 
up to log Luv i= 10.3-10.4 (L Q ) (there is only a "discrepant" 
point at z = 0.7 but with a very large error bar). For the high- 
est observed UV luminosities (log(Luv[L Q ]) > 10.3-10.4), only 
present in the samples at z — 1, the fraction of galaxies with 
log(LiR/Luv) > 0.7 decreases both for the U selection and the 
LBGs. The LBG sample has a slightly lower fraction of galax- 
ies with log(LiR/Luv) > 0.7 than that found in the U selected 
sample over the whole range of luminosity. 

The fraction of galaxies with log(LiR/Luv) > 1 highlights 
galaxies with the largest dust attenuation. This fraction does not 
exceed ^ 20 % for all our samples. Galaxies with such a high 
extinction seem to be more frequent at z > than at z = 0: the 
distribution of Lir/Luv is found to reach larger values at z > 
than at z = but we must remain cautious because on the 
uncertainties on the MIR-total IR calibration. As found above , 
at z = 1 there is almost no UV luminous galaxy with a large 
attenuation and the fraction of galaxies with log(Li R /Luv) > 1 
increases toward lower UV luminosities. 



The evolution of the fraction of LIRGs is reported in Fig. [6 
This fraction increases with the UV luminosity: it is expectec 
even even without any evolution of the Lir/Luv distribution. For 
the UV selected galaxies at z = 0, 0.7 and 1, the variations are 
found similar with a slightly larger fraction of LIRGs for a given 
UV luminosity at z > 0, it is indeed the same effect as noted for 
galaxies with log(Lni/Luv) > 1- The fraction of LIRGs in the 
LBG sample is systematically lower than that found for the UV 
selected galaxies at z = 1, again leading to the conclusion of a 
lower dust attenuation for these galaxies. 

|Reddy et al.| ( 2008] > estimated the colour excess distribu- 
tion of BX galaxies at z — 2 and LBGs at z - 3. They found 
< E(B-V) >= 0.15 + 0.07. Adopting the dust attenuation law of 
|Calzetti et al. (2000) gives < Auv >= 1.5mag. As a consequence 
50% of BX galaxies and LBGs have Auv > 1-5 mag and if we 
assume that the distribution of E(B-V) is Gaussian, 30% have 
Auv > 2 mag. |Reddy et al.| (|2P08 ) obtained similar results by an 
analysis of the L IR /Luv distribution. They also found an aver- 
age dust attenuation which does not vary with the UV rest-frame 
luminosity. 

Therefore dust attenuation in UV selected galaxies at z > 1 
seems to be slightly larger than that found at z — 1, the dif- 
ference being particularly significative for UV luminous galax- 
ies for there is a hint for a lower dust attenuation at z — 1. 
Nevertheless we must remain cautious in our conclusions given 
the uncertainties in the estimates of dust attenuation and the dif- 
ferent methods adopted: IR to UV flux ratio up to z = 1 based on 
12 and 15 /zm luminosities, UV colours and 8 /mi luminosities 
at higher z. These methods are known to give different results at 
least up to z = 1 ( |Burgarella etaL][2007||Elbaz et al. [[20071 ) and 
the calibration of the 8 /mi luminosity in total IR luminosity is 
highly uncertain ( Caputi et al. 2007 Burgarella et al. 2009 [ l. 



3.3. Lir/Luv versus L K 

We have seen that it is difficult to interpret the variation of 
Lir/Luv as a function of a quantity also depending on these two 
luminosities. We can also use a quantity independent of them 
to avoid such an issue. In this section we will consider the rest- 
frame K luminosity of the galaxies which is a tracer of the stellar 
mass of galaxies. They are calculated with the IRAC band at 3.6 
and 4.5 /mi which corresponds to rest-frame K at z — 0.7 and 1 
and with 2MASS data at z = ( |Buat et al.| [2007] l In Fig. [7] is 
reported the variation of Lir/Luv as a function of Lk for the dif- 
ferent samples and also split in luminosity bins. A net increase 
of Lir/Luv with Lk is found without a clear evolution with z 
for the galaxies selected in UV rest-frame, only LBGs appear to 
have a lower dust attenuation for a given Lk . When the sam- 
ples are split according to the UV luminosity of the galaxies it 
appears in average that the more UV luminous objects exhibit a 
lower Lir/Luv for a given K luminosity than UV fainter sources. 
This is in agreement with what has been found in section 3.2. A 
linear regression on both Luv and Lk gives 

log(Li R /L uv ) = 

0.78(0.04) log L K - 0.79(0.06) log L uv + 0.86(0.21) (5) 

Excluding LBGs from the analysis (since they are less extin- 
guished than UV selected galaxies) leads to a slightly different 
regression 

log(L IR /L uv ) = 

0.78(0.04) \ogL K - 0.65(0.07) log L uv - 0.58(0.21) (6) 
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Fig. 7. L IR /L UV versus the K rest-frame luminosity L K expressed in solar units. Left panel: the different samples are plotted with 
the same symbols as in Fig. [2] Right panel: the samples are split according to the UV luminosity as in Fig. [3] 



Using the calibration of Chary & Elbaz (2001 ) would give 



log(L IR /L uv ) = 

0.81(0.04) log L K - 0.54(0.06) log Luv 

for the whole sample and 



1.92(0.26) (7) 



log(Li R /L uv ) = 

0.81(0.04) logLtf - 0.54(0.08) log L uv - 1.91(0.25) 



(8) 



when LBGs are excluded. 

Marti n et al. (2007) and Iglesias-Para mo et al. (2007} also 
studied the variation of L IR /Luv as a function of the stellar mass 
from z — to 1 for UV selected galaxies, our results are globally 
consistent with theirs: the relation between Li R /Luv and Lk or 
M stm - exhibit the same steepness, the shift with redshift reported 
in these studies is of similar amplitude to the one we find as a 
function of Luv- The main differences between these previous 
analyses and ours is that we select galaxies with the same limit 
in UV luminosity whereas the other studies were based on mag- 
nitude limited samples with a different limit in luminosity when 
z varies. Our approach allows us to emphasise differences of be- 
haviour as a function of the luminosity of the galaxies and this 
effect seems to be at the origin of the redshift evolution reported 
earlier. 



4. Star formation activity 

Dust attenuation has been found to be related to the observed UV 
and K luminosity of our UV selected galaxies, LBGs exhibiting 
more extreme behaviour. Is the strength of the star formation ac- 
tivity also linked to the UV luminosity? Do LBGs exhibit differ- 
ent star formation activity than UV selected galaxies? The star 



formation activity can be quantified by the specific star forma- 
tion rate defined as the current star formation rate divided by the 
stellar mass of the galaxies. We can estimate this quantity with 
our data sets. We calculate the total SFR by adding the SFR from 



the IR and the UV (observed) luminosities (Iglesias-Paramo et 
|aL] [20061 |Buat et al.| [20081. We adopt a Sa 
formulae of Iglesias-Paramo et al. ( 2006| l: 



peter IMF and the 



log SFR IR (M Q yr- 1 ) = log L IR (L ) - 9.75 
log SFR uv (M yr- 1 ) = log L UV (L Q ) - 9.51 



(9) 
(10) 



The total SFR is expressed as SFRi R + SFRuv except at z — 
where the contribution of the dust emission not related to the star 



formation is estimated to be 30 % (Iglesias-Paramo et al. 2006 ). 
The stellar masses of the galaxies are calculated with the IRAC 
band at 3.6 and 4.5 fim which corresponds to rest-frame ^T-band 
at z — 0.7 and 1 and with the calibration of Arnouts et al.| ( |2007] l. 
At z = 0, the calibration of Bell et al. (2003 ) is used as discussed 
H2006) and |Buat et al.| (|2008). We adopt 



in 



Iglesias-Paramo et al. 



a Salpeter IMF and check that the extrapolation of the calibration 
of | Arnouts et al.| ( [200 7 ) is consistent with that we adopt at z — 
within 30% (0.1 dex). 

In Fig.[8]is reported the variation of the SSFR as a function of 
the stellar mass for each sample truncated at log Luv > 9.9 (L Q ). 
Such a limit combined with the detection limit adopted for the 
fluxes at 24 fim (25 yt/Jy) gives a limit in SFR for z — and 
1 also indicated in Fig. [8] The SSFR at a given stellar mass 
increases with z as reported in both observational and theoret- 
ical studies and predicted in scenarios of galaxy evolution. The 
SSFR also exhibits a flat distribution: it is expected when only 



star forming galaxies are selected (|Elbaz et al. 2007 Buat et 



[ay [2008 ). The consistency between models and the mean trends 
found in UV and IR selected samples has been shown to be good 
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Fig. 8. Variations of the specific star formation rate for every sample considered in this work, all the samples are truncated at 
log Luv > 9.9 in solar units. Left panel: as a function of the stellar mass of the galaxies . The diagonal lines are the lower limits 
in SSFR for z = 0.7 (lower line) and 1 (upper line). Right panel: as a function of the UV luminosity for galaxies with log(M star ) > 
10.8(M o ) 



up to z = 0.7 ( |Buat et al.| [2008 ) but breaks at z > 1 ( [Elbaz et| 
[ak] |2007| l. The purpose of this paper is not to perform a com- 
parison between models and observations since we are dealing 
with only a sub sample of the overall galaxy population: the ob- 
jects with log Luv > 9.9(L Q ). We want to compare the proper- 
ties of these galaxies at different z and with LBGs. LBGs and 
UV selected galaxies seem to experience similar SSFRs at the 
same redshift. We now focus on the most massive galaxies with 
log(M stal ) > 10.8(M Q ) for which the detection limits reported in 
Fig.[8]do not induce a substantial bias. The galaxies of this sub- 
sample with a moderate UV luminosity (log Luv < 10.4(L Q )) 
exhibit a large range of SSFRs, the most quiescent objects be- 
ing found at z=0. Conversely all the most UV luminous objects 
(log Luv > 1O.4(L )) exhibit quite large SSFRs. These galaxies 
experiment high SFRs between 10 and 190 M yr~' with an aver- 
age value of 55 M yr~' : with such a rate they might have formed 
all their mass in few Gyr. 



5. Total UV+IR luminosity functions 

Do we miss star forming galaxies in a UV selection up to z = 1, 
and as a consequence are we able to retrieve all the star formation 
when applying a reliable dust attenuation to galaxies selected 
in UV rest-frame? To answer these questions, it would be the 
most direct to construct the luminosity functions (LFs) with total 
luminosity related to star formation activity, Luv + Lir. 



5.1. Method 

The most important but difficult point of this analysis is how 
to handle the two variables Luv an d at the same time. We 
proceed a statistical analysis as follows: 

Step 1 Since our sample is primarily selected at UV (GALEX 
FUV at z = taken as reference (A = 1530 A), GALEX 
NUV at z = 0.7 and EIS {/-band at z = 1), we construct 
univariate UV LFs. 

Step 2 We bin the UV LFs, and estimate the distribution of the 
total IR luminosity estimated by equations ^ and (|4]) at each 
bin. 

Step 3 We sum up the distribution functions of the total IR lu- 
minosity along all the UV luminosity bins. 

For Step 1, we used improved versions of two represen- 
tative LF esti mators : l/V max -est imator ( |Schmidt| |1968| l and 
C -estimator (Lynden-Bell 1971 1 in an optimal manner, ex- 



plained and examined extensively by |Takeuchi et al.| ( |2000| i and 
Takeuchi (2000). We estimated the LPs of NUV-selected sam- 
ple at 0.6 <z. < 0.8, {/-band selected sample at 0.8 < z < 1.2, 
and LBG samples 0.9 < z < 1.3. Since we are interested in the 
star-forming galaxies, we have omitted known quasars/AGNs 
from our sample, as already explained in section 2.1. These uni- 
variate UV LFs are estimated simply at their primary selection 
bands for the NUV-selected z = 0.7 sample and {/-band selected 
Z = 1.0 sample: we obtain LFs at ~ 1400 A and ~ 1800 A for 
Z = 0.7 and 1 respectively. As discussed in section 2.1, the NUV 
and U selected samples correspond to rest-frame wavelengths 
close enough to avoid K-corrections. It is not the case for the 
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LBG sample primarily selected in NUV which correspon ds to 
~ 1 100 A rest-frame. We will go back to this issue in Section 5.2 



In Step 2, we should be careful for a significant number 
of upper limits at MIPS 24 /mi. For this, we made use of the 
Kaplan-Meier estimator which enables us to make use of the 
information content earned by the upper limits, originally de- 
veloped in the field of lifetime data analysis (Kapla n~& Meier| 
1958| . Another desirable property of the Kaplan-Meier esti- 
mator is that we can obtain its variance in an analytic form. 
Formulations, derivations and some important properties will 
be discussed and explained elsewhere (Takeuchi et al. 2009, in 
preparation). 

In Step 3, statistical errors are summed up in quadrature, i.e., 
the variance from the primary univariate UV LF and that from 
the IR luminosity distributions. However, we did not include the 
variance caused by the density inhomogeneity of galaxies (often 
referred to as the cosmic variance). We will discuss this issue in 
SectionT 



5.2. UV luminosity functions 

We show the UV LFs for the samples at z = 0.7 and 1 in Figure[9] 
Since both l/V m . dx and C~-estimates agree very well with each 
other, we only show the latter in this paper. We also estimated 
the UV LF of LBGs, but it is not obtained by the same method 
as the other two, since the primary selection is done at NUV 
(for the redshift range of 0.9 < z < 1.3, this corresponds to 
1 100 A). Then, we must use a bivariate method exactly the same 
as described in above Step 2: first we construct a univariate UV 
LF at 1 100 A, and construct distributions of UV luminosities at 
each bin from the {/-band data (~ 1700 A at the rest-frame), and 
sum them up along the UV luminosity bin. The obtained UV LF 
of LBGs is presented in cyan symbols in Figure|9] 

Roughly speaking, we observe that the UV LFs at redshift 
ranges of z — 0.7 and z — 1.0 agree with those of |Arnouts et 
al. (2005). This shows that our selection is appropriate for the 
purpose of this study. We may see that there are some discrep- 
ancies between the shapes of our UV LFs and those of Arnouts 



et al.| (2005 1. The difference may be attributed to the difference 
in selection of galaxies, e.g., Arnouts et al. (2005) performed a 
NUV selection and Zf-corrected the flux, unlike our selection at 
U -band at z — 1 without ^-correction. The different way pho- 
tometry has been performed can also be at the origin of a subtle 
difference: the GALEX deep fields are known to be crowded and 
down to NUV — 23 mag, Sextractor does not separate accurately 
the sources leading to an under-density and a brightening of the 
sources. In the {/-band the PSF is smaller, so less or not affected 
by confusion. 

We then focus on the difference in the UV LF from U -band 
and the LBG LF. As mentioned above, the former was con- 
structed with the [/-band selection with the univariate method, 
while the latter was constructed by the bivariate method: we first 
select galaxies at NUV and with a FUV-NUV colour criterion, 
then we estimate the {/-band luminosity distribution. Therefore 
LBG selection criteria at z = 1 . 1 are not very different from the 
{/-selection at similar redshift (z — 1) since both are based on 
a UV rest-frame selection. Nevertheless, in the LBG selection 
we miss UV-faint galaxies at L uv < 10 10 L Q when compared 
to the {/-selection probably because of the combined effects of 
a selection at a shorter rest-frame wavelength for the LBGs and 
of a FUV-NUV criterion which selects only blue objects, as ex- 
plained in section 2.2. This suggests that the LBG selection cri- 
terion is more restricted to pick up UV-luminous galaxies. We 



will come back to this point when we discuss the difference of 
the UV + IR LF. 



5.3. The total UV+IR luminosity functions 

Here, we show the total UV+IR LFs from our UV-selected sam- 
ples. Again we stress that we addressed the upper limits of the 
sample at MIPS 24 fim by Kaplan-Meier method, i.e., we have 
made a maximal use of the observed information from IR. We 
show the total LFs in Figure 10 Top panel shows the Luv + Lir 
LF at z — 0.7, while bottom panel is the one at z = 1.0. In 



Figure 10 we also show univariate UV LFs constructed from 



purely UV-selected samples by Arnouts et al. ( 2005 1 (dotted 
lines), as well as univariate IR LFs made from purely IR-selected 
samples at 24 //m by Le Floc'h et al. (2005) (dashed lines). The 
symbols are the LFs derived from our sample. Errors are cal- 
culated analytically by the asymptotic variance formula of the 
Kaplan-Meier estimator, convolved with the statistical error of 
the univariate LFs at UV. The indicated errors are lcr (68 % CL). 
Because of the known limitation of the Kaplan-Meier estimator, 
the lowest luminosity bins are underestimated (indicated by ar- 
rows on the symbols). 

Clearly, the total LFs are much higher than the univariate UV 
LFs. This means that most of the luminosity of a galaxy at these 
redshifts is emitted in the IR. Since the luminosity related to their 
star formation activity tends to be emitted in the IR wavelengths 
(e.g., Takeuchi et al. , 2005b ), the resulting total LFs is consistent 
with this known fact. 

At z = 0. 7, the total LF is even higher than the IR LF of|Le] 



Floc'h et al.|p005 



but within a range of the cosmic variance 
(~ 60 % for GOODS: |Somerville eTaLl[2004) l. Apart from this, 
it is rather consistent with the IR LF. This is an expected higher-z 
counterpart of the result discussed by Buat et al. (2007} at z = 
0.0. 

In contrast, the total LF is significantly lower than the IR 
LF at z — 1.0 for galaxies more luminous than 2 x 10 n L o . 
It is worth mentioning that the primary UV LF has an excess 
in normalisation compared with the global univariate UV LF at 
the same redshift. Then, this deficiency of galaxies turns out to 
be quite significant. Though the most luminous bin is disturbed 
by the symbol with a very large error, we see a trend that the 
more luminous galaxies are, the larger the discrepancy becomes. 
This is a clear piece of evidence that our UV-selection misses 
intrinsically luminous galaxies which are active in star forma- 



tion. Such galaxies must be studied from an IR selection. Buat 
et al. ( 2007| l studied both IR and UV selections in the nearby 



universe and found that the more luminous galaxies (in terms of 
total Luv + Lir luminosity) are present in the IR selection and 
suffer a very strong dust attenuation. The galaxies more lumi- 
nous than =i 2 x 10 n L Q exhibit a mean log(LiR/Luv) ~ 2. The 
relation found between Luv + Lir and Lir /Luv at z — has been 



found to be globally valid at higher z for IR selections ( Choi et 



[aL] [2006) |Xu et al.j [20071 |Zheng et al] |2006] >. An L lR /L VY ratio 
larger than 100 for galaxies more luminous than 2 x 10 n L Q 
implies that these galaxies are not detected at z = 1 in our U 
selection (limited to log(Luv[L ]) > 9.9) 

At z - 0.7 our UV selection goes deeper ( log(Luv[L Q ]) > 
9.3), therefore galaxies with a larger dust attenuation can be 
detected in UV. It is also interesting to note that the galaxies 
with the highest Lir /Luv are found at z = 0.7 (Fig. [4} for 
the UV faintest galaxies and that the fraction of galaxies with 
log(LiR/Luv) > 1 is globally higher at z — 0.7 than at any other 
redshift (Fig. [5]). However, it is puzzling to invoke a drastic evo- 
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lution from z = 0.7 to z — 1 .0. Since the cosmic time differs by 
less than a few Gyr, such an evolution should be very fast. 

As a conclusion, up to z = 1 UV rest-frame observations 
must be much deeper (by more than 5 mag) than the expected 
limit in bolometric luminosities if one must be able to detect 
most of the star forming galaxies. 



Figure 11 shows the total UV+IR LF of the LBG sample. 



The deficiency of total LF is more prominently seen in the 
LBG LF. In this case, the deficiency ranges toward lower lu- 
minosities a 4 x 10 I0 L o . Considering the LBG sample selection 
which makes use of NUV and FUV fluxes observed at 2310 and 
1530 A by GALEX, this trend may be understood consistently: 
the LBG sample consists of galaxies with less extinction on av- 
erage, leading to less IR luminosities with respect to the same 
Luv as shown in section 3. The deficiency with respect to the 
{/-selection affects the faintest bins of the LBG LF (Fig. [9] and 
discussion in section 5.2). Since the dispersion of the Lir/Luv 
distribution is very large as we have seen before, the contribu- 
tion of these bins to the number density of galaxies is significant 
and the deficiency of galaxies affects all the range of the total 
UV + IRLF. 



5.4. Discussion 



As discussed by many authors ( e.g., |Chary & Elbazj [2 001 ; 
|Takeuchi et aTl|2005a||Caputi et al.[|2007||Rieke et al.||2008| ), the 
monochromatic MIR luminosity-to-total IR luminosity conver- 
sion plays an important role, especially because of the currently 
limited number of deep multi-band observations at FIR. Since 
the intrinsic scatter in the linear regression is not very small, 
these "calibration formulae" inevitably have significant uncer- 
tainty. Then, it makes sense to examine how different formulae 
change our results, especially the deficiency of intrinsically lu- 
minous galaxies in the UV-selection. 

To test this, we have estimated the total UV + IR LFs exactly 
in the same manner but with CE01 conversion. The resulting LFs 
are shown in Figure 12 In Figure 12 open squares represent the 
LFs with the formula of Takeuchi et al. (2005a), while open tri- 
angles are the ones with CE01 conversion. All the other symbols 
are the same as in Figure [TO] 

At z = 0.7, since the difference of these formulae is quite 
small (cf. upper panel in Fig. |TJ, the results are almost the same. 
At z = 1.0, the difference is visible between the two estimates. 
As we have seen in Figure [T] CE01 formula gives larger IR lu- 
minosity. Hence, it produces larger total luminosity in Luv + Lir 
with respect to the same L uv . As a result, the discrepancy be- 
tween the total LF and the IR LF becomes smaller, but still sta- 
tistically significant. 

At z - 2 and 3, |Reddy et al. ( |2008[ ) were able to reproduce 



all the IR LF up to L IR = 10^ L Q from only UV-optical data with 
even an excess of faint sources as compared to the results from 
IR surveys alone. They built the UV rest-frame LF with Monte 
Carlo simulations to recover all the star forming galaxies; then, 
to recover the IR LF they assumed either a constant dust attenu- 
ation distribution irrespective of UV luminosity or a decrease of 
the average dust attenuation for UV faint galaxies. Conversely, 
at z - 1 dust attenuation is found not very dependent on UV 
luminosity with only a slight decrease for UV luminous galax- 
ies. We cannot reconstruct the bright end of the IR LF from a 
UV-selected sample. The IR luminous galaxies, observed in IR 
surveys, exhibit a very large dust attenuation which makes them 
undetected in UV (rest-frame). Since we use actually observed 
IR and UV flux densities (including upper limits) our method can 
be considered as being secure, although it is dependent on the va- 



lidity of the MIR to total IR luminosity conversion. The method 
of Reddy et al. also suffers from the uncertainty in the MIR to to- 
tal IR luminosity conversion which is particularly large at 8 /im 
rest-frame and on the accuracy of dust attenuation factors esti- 
mated from the UV-optical alone (see discussion in section 3.2). 
Nevertheless if we trust both results it implies a lower fraction 
of galaxies intrinsically UV+IR luminous and with a large dust 
attenuation at z — 2-3 than at z - 1. We will re-investigate 
this issue by using IR-selected samples up to z = 1 in a fully bi- 
variate manner (Takeuchi et al. 2009, in preparation). The future 
observations of HERSCHEL should give us the high redshift IR 
selected samples necessary to solve this question. 



6. Conclusions 

We have analysed the IR emission of galaxies selected in UV 
rest-frame from z = to z = 1. The samples were built to be 
very homogeneously selected in terms of wavelength and lumi- 
nosities. We also considered a sample of Lyman Break Galaxies 
atz =i 1. 

1. The Lir/Luv ratio was used as a proxy for dust attenua- 
tion. For the bulk of our galaxy samples, this dust attenua- 
tion is found not to vary significantly with z and less than 
20 % of the sample galaxies have a Lir/Luv > 10. Anyway, 
some evolution is seen in the extreme regimes of high and 
low Lir/Luv ratio. The most luminous UV objects (Luv — 
2 x 10 10 L Q ) are only present at z — 1 and exhibit a very low 
dust attenuation. When Luv S= 2 x 10 10 L o , the fraction of 
galaxies with a high Lir/Luv is larger at z > than in the 
nearby universe and the galaxies with the largest dust atten- 
uation are the faintest ones in our samples (Luv — 3 X 10 9 L Q ) 
. However these results all depend on the MIR-total IR cal- 
ibration which is uncertain. Dust attenuation increases with 
the K luminosity in a similar way at all redshifts. A residual 
trend is found with the UV luminosity: when Luv increases, 
galaxies of a given Lk have a lower Lir/Luv- A relation be- 
tween Lir/Luv, Luv an d L# is given. LBGs at z — 1 seem to 
be less extinguished than UV selected galaxies of similar UV 
luminosity and at same z. Since the UV luminosity of galax- 
ies globally increases with z, these trends found with the UV 
luminosity must be accounted for to interpret the evolution 
with redshift of Lir/Luv reported in previous studies. 

2. Massive and UV luminous galaxies (log(M star ) > 10.8(M Q ) 
and log Luv > 1O.4(L )) are found very active in star forma- 
tion (large SSFR) whereas fainter galaxies of similar mass 
exhibit a larger range of SSFR. LBGs and UV selected galax- 
ies have similar SSFR. 

3. We constructed LFs with total luminosity related to star for- 
mation activity, Luv + Lir from our UV-selected galaxy sam- 
ples. We have used the Kaplan-Meier estimator to make use 
of information carried by IR detections and upper limits in a 
coherent manner. The resulting total UV + IR LFs are much 
higher than the univariate UV LFs from purely UV-selected 
samples. This means that most of the luminosity produced 
by star formation activity is emitted in the IR wavelength 
range. Though at z — 0.7, the total LF we obtain is consistent 
(even higher because of a density excess) with the univariate 
IR LF, we find a clear deficiency of galaxies in the total LF 
at z = 1.0 for galaxies more luminous than 2 x 10 n L Q 
This result is not significantly affected by different total IR 
luminosity calibration formula. Thus, we conclude that the 
IR LF cannot be reconstructed solely from our UV-selected 
galaxies at z = 1 and that deeper data are needed in order 
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to detect galaxies with a large Lir/Luv- Practically, to detect 
most of the star forming galaxies down to a given bolomet- 
ric magnitude, UV rest-frame observations must be deeper 
than this bolometric limit by at least 5 mag (corresponding 
to (Lir/Luv — 100). The deficiency in the total LF is found 
much higher for the LBG selection affecting the whole range 
of luminosity explored in this work for these objects (i.e. 
> 4 x 10 I() L Q ). 
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